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I.  INTRODUCTION 


Accurate  modeling  of  the  HF  and  DF  chemical  lasers  requires  a 

reliable  set  of  reaction  rate  coefficients.  Cohen  *  has  reviewed  the  literature 

to  provide  a  standardized  rate  package  for  use  with  computer  codes.  In 

addition  an  excellent  review  of  gaseous  fluorine  reactions  has  been  recently 

2 

presented  by  Foon  and  Kaufman.  The  rates  of  the  so-called  "cold"  reactions, 
F  +  and  F  +  D2«  are  among  the  most  important  parameters  in  the  laser 
systems: 


f+h2— Lhf+  +  h  (1) 

F  +  D2— 5»DFt  +  D  (2) 

A  number  of  measurements,  both  relative  and  absolute,  have  been  reported 

3 

in  the  last  several  years.  Just  recently.  Quick  and  Wittig  and  Preses 

4 

et  al.  reported  that  the  infrared  multiphoton  dissociation  of  SF^  is  a  con¬ 
venient  F -atom  source  for  time -resolved  kinetic  measurements.  In  the 
present  work  this  technique  was  used  under  demonstrably  isothermal  experi¬ 
mental  conditions  to  obtain  extremely  precise  rate  coefficients  for  these  two 
reactions.  The  accuracy  of  the  determinations  critically  depends  on  the  use 
of  the  total  HF(DF)  fundamental  chemiluminescence  as  a  diagnostic  for  the 
extent  of  reaction.  This  assumption  is  examined  both  analytically  and  with 

5 

the  use  of  The  Aerospace  Corporation  numerical  modeling  code  NEST. 
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H.  EXPERIMENTAL  APPARATUS  AND  PROCEDURE 


The  measurements  were  performed  in  a  room -temperature  flow 
system.  Argon  and  calibrated  mixtures  of  SF^  in  argon  (10.24%),  in 
argon  (10.  14%),  and  in  argon  (10.  37%)  were  of  the  highest  purity  available 
in  bulk-gas  form.  Molar  flow  rates  were  measured  with  rotameter  flow¬ 
meters  calibrated  against  pres  sure- rise  measurements  in  a  standard  volume. 
The  total  pressure  was  measured  with  a  capacitance  manometer  gauge 

_3 

(MKS  Tru-Torr)  with  a  resolution  of  1  X  10  Torr.  All  data  were  taken  at 

_3 

3.  95  ±  0.  03  Torr  pressure.  SF^  partial  pressures  of  6.  5  X  10  and 

_3 

33  X  10  Torr  were  used,  and  the  H2CD2)  partial  pressure  was  varied  from 

.3 

2.5  to  75  X  10  Torr.  Argon  constituted  the  balance  of  the  flow. 

The  fluorescence  cell  incorporated  into  the  flow  system  was  a 

10-cm-i.d.  brass  cube  internally  coated  with  teflon  (Fig.  1).  Radiation 

from  a  pulsed  CC>2  laser  (Lumonics  K-202-2)  entered  and  exited  the  cell 

through  15-cm-long,  2.5-cm-i.d.  side  arms  fitted  with  NaCl  windows .  The 

laser  was  focused  with  a  38-cm  focal  length  f .  1.  ZnSe  lens  into  the  center  of 

the  cell.  In  this  region  the  beam  diameter  was  constant  at  approximately  2  mm 

over  5  cm  of  optical  path  length.  The  reported  data  were  taken  with  the  use 

of  the  P(20)  CO^  line  at  10.  59  pm  at  an  energy  of  1. 5  J /pulse.  No  visible 

emission  was  detected  with  a  GaAs  photomultiplier  when  argon  alone* 

Ar  +  SF^,  or  Ar  +  was  flowed  through  the  cell.  Thus  gas  breakdown  was 

eliminated  as  a  consideration  in  these  experiments.  Substantial  HF  3  -*  0 

o 

overtone  emission  was  observed  at  8900  A  when  and  SF^  were  simul¬ 
taneously  present. 


SIGNAL 

AVERAGER 


TRANSIENT 

RECORDER 


+  + 

The  chemiluminescence  from  HF  or  DF  was  collected  by  an  internal 
2-in.  -diameter  i/2  aluminum  mirror  overcoated  with  MgF£.  The  radiation, 
after  appropriate  filtering,  was  focused  onto  a  77  K  InSb  detector  (Texas 
Instruments).  The  amplified  signal  was  recorded  with  a  Biomation  805 
transient  recorder  and  stored  in  a  Nicolet  1072  digital  signal  averager. 

The  data  represent  an  average  of  eight  laser  pulses.  A  typical  trace  is 
shown  in  Fig.  2  with  the  rise  and  fall  of  the  chemiluminescence  recorded  on 
separate  time  bases.  The  measured  rise  time  was  always  greater  than  10 
times  the  detector -amplifier  response  time  of  1.  5  psec.  The  gas  mixture 
was  replaced  approximately  five  times  between  laser  shots  in  order  to 
eliminate  the  possible  effects  of  product  accumulation  on  the  measured  rates 
Although  the  experimental  arrangement  is  nominally  a  flow  system,  the 
measured  fluorescence  times  are  orders  of  magnitude  shorter  than  the 
transit  time  of  gas  through  the  cell.  Thus  the  data  are  analyzed  in  terms 
of  a  laser -irradiated  static  mixture. 
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.HR**-, 


t  (arbitrary  units) 


Fig-  2.  Time -Resolved  Infrared  Chemiluminescence  Trace. 

Channels  1-256,  0.5  (as ec/ channel;  Channels  513-1024, 
10  (jisec /channel.  Partial  Pressures:  SF&  =  33  mTorr, 
=  24.3  mTorr,  Ar  =  3.87  Torr. 
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III.  DATA  ANALYSIS 


The  rate  coefficient  kj  can  be  determined  directly  from  the  rise  time 
of  the  fluorescence  from  HF^  produced  by  Reaction  (1).  This  technique  has 
been  used  previously  for  the  measurement  of  reaction  rate  coefficients  for 
reactions  of  Cl  atoms  with  HBr(DBr)  and  HI(DI).  ^  It  provides  an  accurate 
value  of  the  rate  coefficient,  if  certain  spectroscopic  and  kinetic  conditions 
are  met.  These  conditions  are  described  in  this  section.  The  total  fluores¬ 
cence  I(t)  is  related  to  the  individual  vibrational  level  populations  by 

I(t)  =  E  Av[HF(v)]  (3) 

th 

where  A  is  the  radiative  lifetime  of  the  v  vibrational  level.  The  initial 
v 

vibrational  populations  are  produced  by  Reaction  (1): 

k.(v) 

F  +  H2  HF  (v)  +  H,  v  =  1,  2,  3 

These  levels  are  removed  by  a  variety  of  processes;  however,  the  following 
analysis  focuses  on  single -quantum  deactivation  by  quenching  partners  such 
as  and  argon,  whose  concentration  remains  essentially  constant  during 
the  experiment  ([F]«  [h2]). 

k4(v) 

HF(v)  +  M  HF  (v  -  1)  +  M  (4) 

Process  (4)  can  populate  certain  of  the  lower  vibrational  levels. 
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Radiative  loss  of  HF^,  V  -  V  coupling  between  excited  HF  levels,  and 
back  transfer  of  vibrational  energy  from  H2(v)  to  HF (v' )  are  omitted  from 
the  analysis.  The  result  is  an  analytic  expression  that  is  tractable.  The 
omitted  terms  were  included  in  the  NEST  computer  calculations  and  found 
to  have  little  quantitative  effect  for  the  experimental  conditions  of  this  study. 
In  this  context,  the  rate  equation  for  [HF (v)]  is 

gf[HF(v)]  =  kt(v)[H2][F]  -  k4(v)[HF(v)][M]  +  k4(v  +  l)[HF(v  +  1 )]  [M] 

(5) 

and  the  rate  equation  for  total  fluorescence  intensity  can  be  written  as 

=  ^2Av[HF(v)]  =  [H2][F]Lk1(v)Av  -  Lk4(v)Av[HF(v)][M] 

+  Lk4(v  +  l)Av[HF(v  +  I)]  [M]  (6) 

Equation  (6)  reduces  exactly  to  Eq.  (7) 

=  [H2][F]LAvk1(v)  -  k4(l)[M]l(t)  (7) 

when  the  deactivation  rate  coefficients  and  the  Einstein  coefficients  are 
related  by  Eq.  (8). 

k4<2>  A2 

k4(l)  A2  -  A j 

k4(3)  A3 

kJT)  =  a3  -  a2  (8) 
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In  the  case  of  harmonic  oscillators,  Av  <x  v  and  kv  oc  v  such  that  Eq.  (8)  is 
satisfied  exactly.  HF  is  not  harmonic,  so  k,j(v)  must  increase  as 
i/2.4/4.93  to  satisfy  Eq.  (8)  with  accepted  values  of  the  Einstein  coefficients. 
Since  [F]  =  [F]q  exp( -k^ [H^]^ ),  Eq.  (7)  can  be  integrated  to  yield 

Ck  [H2][F] 

I(t)  =  k1[H2]-k4(l)[M]  Je*P(-k4(l)[M]t)  -  expf-kjHLjt)'  (9) 


where  C  is  given  by  (LA^kj  (v))/k^ .  Therefore,  the  rise  time  of  the 
fluorescence  intensity  is  directly  related  to  the  rate  coefficient  by 
=  (k 4 [H^ ] )  *.  snd  the  decay  time  can  be  expressed  as  =  (k4(l)[M])~*. 
Equation  (9)  can  be  rewritten  as; 


I(t)  = 


1 


.-l 


exp 


(10) 


At  long  times  the  intensity  is  represented  by 


iL(ti  = 


cT;‘tF)o 


"r  Lrj0  /-t\ 


(ID 


The  short-time  intensity  difference  between  Eqs.  (10)  and  (11)  is  given  by 


Ct_1  [F]  /  v 

Is(t)  =  IL(t)  -  !(t)  =  - rfexp  “ 

T  -  T,  \  Tf 


(12) 


Both  Eq.  (11)  and  Eq.  (12)  have  the  same  intercept,  which  can  be  related  to 
the  initial  F-atom  concentration 
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[FJn- 


■l'0*  t;‘  -  t;‘ 


o'  c 


is(°)  t;1  -  r-‘ 


(13b) 


Numerical  modeling  calculations  were  performed  with  the  NEST3 
computer  code  with  input  parameters  corresponding  to  the  present  experi¬ 
mental  conditions  in  order  to  test  the  accuracy  of  these  analytic  approxima¬ 
tions.  Recommended  values  of  A  ,  the  reaction  rate  coefficients  k,(v)  for 
the  F  +  H 2  ■•HF(v)  +H  reactions,  and  the  HF(v)-HF  deactivation  processes 
IC4 (v)  were  taken  from  Reference  1.  Rate  coefficients  for  HF(v)  deactivation 
by  H  and  H 2  were  taken  from  References  7  and  8,  respectively. 

The  computer  code  calculated  a  time -dependent  fluorescence  profile 
from  which  a  rise  time  and  a  decay  time  were  extracted  by  the  same  method 
applied  to  the  laboratory  data.  The  limitations  of  the  analytical  equations 
were  revealed  by  comparing  the  derived  rate  coefficients  with  the  input  values 
The  most  important  conclusion  is  that  =  (kjfp^])  *  is  a  good  approximation 
even  when  Eq.  (8)  is  not  valid,  as  long  as  «  t^.  The  computed  fluore¬ 
scence  prof  les  obey  this  latter  criterion  to  a  high  degree  of  accuracy. 

Results  of  these  calculations  are  shown  in  Table  I  for  three  different 
pressures. 


The  experimental  decay  times  are  somewhat  shorter  than  those 


from  6  to  30  for  the  experimental  data.  The  measured  decay  times  were 
sensitive  to  gas  purity  and  handling  techniques;  impurity  deactivation  of 
HF(v),  e.g.,  by  1^0,  may  play  a  role  in  these  experiments.  Artificially 
fast  decay  rates  were  inserted  into  the  calculations  in  order  to  reproduce 
the  experimental  decay  times  and  determine  the  effect  of  additional  deacti¬ 
vation  on  the  derived  rise  times.  Three  methods  of  scaling  this  increased 

2 

deactivation  with  v  were  examined:  (1)  k^(v)  =  constant,  (2)  k^(v)av  ,  and 
(3)  k4(v)  scaling  by  Eq.  (8).  The  normalized  rise  times  and  decay  times 
determined  from  the  fluorescence  profiles  are  listed  in  Table  I  for  these 
three  cases.  The  first  two  cases  did  not  decay  with  single  exponential  rates; 
therefore  the  decay  times  are  approximate  force  fits.  The  calculated  rise 
time  for  0.001  Torr  F  and  0.005  Torr  was  lengthened  somewhat  by 
the  removal  of  20%  of  the  during  the  chemical  reaction.  From  these 
calculations  the  systematic  error  in  relating  the  fluorescence  rise  time  to 
the  rate  coefficients  k^  and  k^  is  estimated  to  be  ±5%. 


Table  I.  Numerical  Modeling  Results  for  Representative 
Conditions  of  0.  001  Torr  F  and  4.  0  Torr  Ar  + 


Case 

0.  005  Torr  H2 

0.  025  Torr  H2 

0.  10  Torr  H2 

la 

1. 00 

0.99 

0.98 

T;'/k,[H2i 

2b 

3 

0.  85 

1. 09 

0.  97 

1.  12 

0.98 

1. 06 

4 

0.  98 

0.  96 

0.95 

1 

na 

na 

na 

1/k4(l)[Ar] 

2 

1.00 

1.  00 

0.98 

I  TC 

3 

1. 02 

1.  14 

1.05 

4 

0.  69 

0.  69 

0.60 

aCase  1  =  Standard  reaction  rate  coefficients  (see  text). 


Case  2  =  k^(v)  scaled  from  Eq.  (8).  1/2.4/4.93/14.3. 

Case  3  =  k4(v)  a  v2. 

Case  4  =  k^(v)  =  constant. 

**In  Cases  2  through  4,  k^(l)  is  adjusted  to  provide  reasonably  close 
agreement  between  calculated  and  measured  decay  times. 


IV.  RESULTS 


The  analytic  description  of  the  infrared  chemiluminescence  is 
approximate,  and  numerous  checks  on  the  accuracy  of  the  analysis  were 
performed.  The  linearity  of  lnl^(t)  and  lnlg(t)  versus  t  is  shown  in  Fig.  3, 
which  represents  a  deconvolution  of  the  trace  presented  in  Fig.  2.  The 
deviations  from  linearity  predicted  in  Table  I  for  certain  deactivation 
rate  scalings  are  barely  discernible  within  the  experimental  scatter.  The 
intercepts  1^(0)  and  lg(0)  were  typically  the  same  (within  10%),  as  predicted. 
Plots  of  [F]q,  calculated  with  Eq.  (13),  versus  [^J  or  ]  were  constant 
to  ±15%  for  nominally  constant  laser  conditions.  Thus  SF^  +  H-^tD^)  colli- 
sions  do  not  contribute  to  [F]q  or,  by  implication,  to  the  time  history  of  HF  . 
These  results  agree  with  previous  deductions  about  the  system.  ^  Although 
each  data  set  was  collected  for  constant  [SF^],  a  plot  of  ln[F]g  versus 
ln[SF^]  had  a  slope  of  1.  06,  demonstrating  that  the  formation  of  F  atoms 
is  linearly  dependent  on  the  [SF^].  As  expected,  SF^-SF^  collisions  do 
not  play  a  role  in  F-atom  formation. 

The  fall  time  of  the  infrared  chemiluminescence  was  plotted  versus 
the  [H^]  and  [D^]  and  compared  with  computer  simulations  that  were  based  on 

*4* 

the  known  collisional,  radiative,  diffusional,  and  convective  losses  of  HF 
and  DF^  (Fig.  4).  The  last  three  removal  processes  are  held  constant  in 
these  plots.  Thus  the  slope  should  represent  the  collisional  quenching  of 
the  coupled  HF(v)  levels  by  H2  and  the  DF(v)  levels  by  D^.  As  discussed 
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Fig.  3. 


Plots  of  1^  Versus  t  [(Eq. 
from  Trace  in  Fig.  2 
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earlier,  the  measured  quenching  rates  are  faster  than  the  computed  rates. 

A  low-level  H^O  impurity  in  the  H^D^)  is  the  most  likely  explanation  for 
these  results,  but  additional  experiments  are  in  progress. 

The  inverse  rise  time  *  of  the  infrared  chemiluminescence  is  given 
by  and  ^[02].  These  quantities  are  plotted  versus  and  [D2I  in 

Fig.  5.  An  unweighted  linear  least -squares  fit  to  these  data  results  in  the 
rate  coefficients  kj  =  (1.81  ±  0.07)*  10*  ^  and  k2  =  (9.  46  ±  0.  30)  X  10*^  cm^/ 
mol -sec.  The  values  given  in  Tables  II  and  III  have  stated  errors  that  reflect 
the  statistical  error  previously  given  and  a  systematic  error  of  ±15%. 

The  determination  of  kj/k2  (Table  III)  is  extremely  insensitive  to 
systematic  errors  in  the  various  calibration  procedures.  The  proposed 
value  kj/k2  =  1.  92  ±  0.  23  has  an  error  given  by  (o  ^  /kj  +  °2^k2  +  05), 

where  a  j  and  o ^  are  the  standard  deviations  of  kj  and  k 2  (Fig.  5).  The 
term  0.  05  ^^2  reflects  the  uncertainty  in  the  use  of  infrared  chemilumine¬ 
scence  from  HF  or  DF  as  an  F-atom  disappearance  diagnostic. 
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Inverse  Rise  Times  tJ  for  (a)  HF  as  Function 
of  Hg  Pressure  and  (b)  DF+  as  Function  of  D ^ 
Pressure 


Table  II.  Absolute  Rate  Coefficients  for  F  +  H,  -»  HF  +  H  at  T  =  298  K 


Table  III.  Absolute  and  Relative  Rate  Coefficients  for 


Multiphoton  Dissociation  (SF 


V.  DISCUSSION 


The  result  for  kj  obtained  in  this  study  is  in  excellent  agreement  with 

3  4 

those  determined  by  Quick  and  Wittig  and  by  Preses  et  al.  Although  in 
all  three  studies  the  same  basic  experimental  technique  was  used,  in  the 
present  study  an  attempt  was  made  to  establish  well-defined  conditions  for 
kinetic  rate  measurements.  Unlike  the  earlier  studies  data  were  collected 

for  a  constant  pressure  of  inert  buffer  gas  (argon)  with  a  ratio  of  [At]/[F]q 

3  24 

~4  X  10  ,  which  eliminates  the  possibility  of  hot -atom  reactions  between 

F  and  H^ID^)  and  limits  the  adiabatic  temperature  rise  to  approximately 

2.  5°C.  Computer  calculations  made  with  the  use  of  the  NEST  code  indicate 

that  the  actual  temperature  rise  during  the  fluorescence  rise  time  is<l  C  for 

[F]Q  =  0.001  Torr.  The  change  in  k^  when  [SF^]q  was  increased  from  6.6  to 

33  mTorr  was  not  statistically  significant.  This  observation  rules  out  a 

temperature  rise  produced  by  the  laser  heating  of  nondissociating  SF^. 

The  other  experimental  determinations  of  k^  (Table  II)  have  been 

1  2 

critically  reviewed  by  Cohen  and  by  Foon  and  Kaufman.  Each  review 

13  3 

recommends  the  value  kj  =  (1.  5  ±  0.  5)  X  10  cm  /mol-sec  and  cites  the 

14 

study  by  Clyne  et  al.  as  the  most  direct.  The  present  measurement  falls 

within  the  high  side  of  that  experimental  limit  and  disagrees  principally  with 

13  15 

several  precise  ESR  studies. 

Absolute  measurements  of  k£  were  nonexistent  until  very  recently 

(Table  III);  thus  the  principal  comparison  of  the  multiphoton  result  must  be 

2 1 

with  earlier  measurements  of  the  kj/k£  ratio.  The  study  by  Persky  is 


the  most  precise  determination  of  kj/l^  made  under  demonstrably 

thermalized  conditions.  That  study  is  in  excellent  agreement  with  the 

23 

nuclear- recoil  experiments  of  Grant  and  Root  in  which  highly  moderated 
1 8 

F  atoms  were  used.  The  measurements  obtained  in  this  study  are  in  close 

agreement  with  these  results  and  are  within  the  experimental  error  of  most  of 

the  other  determinations  in  Table  III.  Thus  the  value  of  =  (9.  5  ±  1.7) 

12  3 

X  10  cm  /mol-sec  is  both  a  direct  determination  and  one  that  satisfies  well- 
established  values  for  k^/k^. 

2 1 

The  temperature  dependence  of  kj/k£  has  been  measured  by  Per  sky, 
[kj/k^  =  (1. 04  ±  0.  02)exp{370  ±  10)/RT],  and  by  Grant  and  Root,  ^ 

[kj/k^  =  (1.04  ±  0.  06)exp(382  ±  35)/RT],  with  upper  temperature  limits  of 
417  and  47  5  K,  respectively.  These  data  provide  a  useful  standard  for 
evaluating  the  results  of  future  high-temperature  studies  of  F  +  and 
F  +  by  means  of  the  multiphoton  dissociation -infrared  chemiluminescence 
technique. 


VI.  CONCLUSIONS 


Accurate  and  precise  rate  coefficients  for  the  F  +  and  F  + 
reactions  have  been  measured  under  highly  controlled  experimental  condi¬ 
tions.  The  extension  of  the  work  to  temperature -dependent  studies  over  the 
wide  temperature  range  required  for  modeling  of  electron -beam -initiated 
HF  and  DF  chemical  lasers  is  in  progress. 
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LABORATORY  OPERATIONS 


The  Laboratory  Operatlona  of  The  Aerospace  Corporation  la  conducting 
experimental  and  theoretical  Investigations  necessary  for  the  evaluation  and 
application  of  scientific  advances  to  new  ellltary  concepts  and  ayeteae.  Ver¬ 
satility  and  flexibility  have  been  developed  to  a  high  degree  by  the  laborato¬ 
ry  personnel  in  dealing  with  the  aany  probleaa  encountered  in  the  Nation's 
rapidly  developing  space  syataas-  Expertise  in  the  latest  scientific  develop¬ 
ments  la  vital  to  the  accomplishment  of  tasks  related  to  these  problems.  The 
laboratories  that  contribute  to  this  research  are: 

Aerophyalca  Laboratory:  Aerodynamics;  fluid  dynamics;  plasnadynanlcs; 
chemical  kinetics;  engineering  mechanics;  flight  dynamics;  heat  transfer; 
high- power  gas  lasers,  continuous  and  pulsed,  II,  visible,  UV;  laser  physics; 
laser  resonator  optics;  laser  effects  and  countermeasures. 

Chemistry  and  Physics  Laboratory:  Atmospheric  reactions  and  optical  back¬ 
grounds;  radiative  transfer  and  atmospheric  transmission;  thermal  and  etate- 
a pacific  reaction  rates  in  rocket  plumes;  chemical  thermodynamics  and  propul¬ 
sion  chemistry;  laser  Isotope  separation;  chemistry  and  physics  of  particles; 
space  environmental  and  contamination  effects  on  spacecraft  materials;  lubrica¬ 
tion;  surface  chemistry  of  Insulators  and  conductors;  cathode  materials;  sen¬ 
sor  materials  and  sensor  optics;  applied  laeer  spectroscopy;  atomic  frequency 
standards;  pollution  and  toxic  materials  monitoring. 

Electronics  Research  Laboratory:  Electromagnetic  theory  and  propagation 
phenomena;  microwave  and  semiconductor  devices  and  Integrated  circuits;  quan¬ 
ta  electronics,  lasers,  and  electro-optics;  communication  sciences,  applied 
electronics,  superconducting  and  electronic  device  physics;  millimeter-wave 
and  far- Inf rated  technology. 

Materials  Sciences  Laboratory:  Development  of  new  materials;  composite 
materials;  graphite  and  ceramics;  polymeric  materials;  weapons  effects  and 
hardened  materials;  materials  for  electronic  devices;  dimensionally  stable 
materials;  chemical  and  structural  analyses;  atresa  corrosion;  fatigue  of 
metals. 

Space  Sciences  Laboratory:  Atmospheric  and  Ionospheric  phyelcs,  radia¬ 
tion  from  the  atmosphere,  density  and  composition  of  the  atmosphere,  aurorae 
and  alrglow;  magnetoapherlc  physics,  cosmic  rays,  generation  and  propagation 
of  plasma  waves  In  the  magnetosphere;  solar  physlca,  x-ray  astronomy;  the  effects 
of  nuclear  explosions,  angnetlc  storms,  and  solar  activity  on  the  earth's 
atmosphere,  Ionosphere,  and  magnetosphere;  the  effects  of  optical,  electromag¬ 
netic,  and  particulate  radiations  In  apace  on  space  systems. 


